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Random optical fields with two widely different correlation lengths generate far field speckle spots 
that are themselves highly speckled. We call such patterns speckled speckle, and study their critical 
points (singularities and stationary points) using analytical theory and computer simulations. We 
find anomalous spatial arrangements of the critical points and orders of magnitude anomalies in 
their relative number densities, and in the densities of the associated zero crossings. 



I. INTRODUCTION 

The canonical exemplars of highly structured light 
fields (complex light) are the random speckle patterns 
produced by scattering of coherent light from surface or 
volume diffusers. Such patterns are not only of spe- 
cial interest in their own right - in part because anything 
generic that can happen in an optical field is likely to 
happen in a speckle pattern - but also because of the 
myriad areas of practical application in which these pat- 
terns appear [1]. 

Prior theoretical and experimental studies of singu- 
larities in speckle patterns have concentrated on fields 
with a single characteristic length scale [1]. Here we 
study speckle patterns with two widely different length 
scales. We call such patterns "speckled speckle" , because 
as shown in Fig. [T] the major speckle spots in the pat- 
tern are themselves highly speckled. We study the basic 
statistical properties of the singularities of these fields, 
in particular their relative number densities, and associ- 
ated level crossings, and find anomalous, often surprising 
results. 

The plan of this paper, which extends a previous brief 
report [2], is as follows. In Section II we describe the 
composite source distributions that generate the nine dif- 
ferent speckled speckle fields studied here, and calculate 
the autocorrelation functions of these fields. In Sec- 
tion III we discuss phase vortices and phase extrema in 
scalar (one polarization component) paraxial fields, and 
C points and azimuthal extrema in vector (two polar- 
ization component) fields. In both cases we find highly 
anomalous ratios for the number densities of singularities 
and extrema - ratios that can differ by many orders of 
magnitude from those in normal speckle. In this section 
we also consider level crossings of the derivatives associ- 
ated with these phase critical points, again finding large, 
orders of magnitude, anomalies for the relative number 
densities. In Section IV we discuss extrema, umbilic 
points, and the associated level crossings, of the real and 
imaginary parts of the field, finding, once again, anoma- 
lously large ratios for these number densities. We briefiy 
summarize our main findings in concluding section V. 

Throughout, wc use computer simulations to provide 
images of the underlying field structures. Based on these 
simulations, which use realistic source parameters, we 
conclude that the many unusual properties of speckled 



speckle described here theoretically could also be studied 
experimentally. 




(b) 

FIG. 1: Speckled speckle, (a) Scalar field (single polariza- 
tion component) intensity, (b) Vector field (two polarization 
component) intensity. The computer simulated source distri- 
butions that produce these calculated speckle patterns consist 
of two concentric disks of light that illuminate a random phase 
plate (diffuser). The diameters of these disks differ by a fac- 
tor of ten, and the total optical power in the smaller disk is 
ten times the power in the larger disk. For the scalar field in 
(a) the illuminating disks have the same linear polarizations, 
for the vector field in (b) the polarizations are orthogonal. 
Individually, the small (large) diameter disk produces large 
(small) speckle spots. Collectively, in (a) interference between 
the speckle fields generated by the small and large disks pro- 
duces the fine granulation of the large speckle spots, whereas 
the pattern in (b) is the sum of the speckle intensities gen- 
erated by each disk because orthogonal polarizations do not 
interfere. This absence of interference notwithstanding, as 
discussed in Section III, the polarization singularity structure 
of the vector field in (b) differs fundamentally from the corre- 
sponding singularity structure produced by each disk acting 
alone. 
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II. SOURCES AND AUTOCORRELATION 
FUNCTIONS 

We assume the (standard two dimensional) geometry 
shown in Fig. [2] further assume circular Gaussian statis- 
tics [1] for the two dimensional, paraxial, complex opti- 
cal field E (x, y) and its derivatives, and write E (x, y) in 
terms of its real TZ (x, y) and imaginary X (x, y) parts as 

E{x,y)^n{x,y)+il{x,y). (1) 

Invoking stationarity and circular symmetry, the normal- 
ized autocorrelation function of the field in the plane of 
observation (the xy-plane) is defined by 



where e in Eq. (6c) is explained below. Here and 



W{v) = {E* {Q)E{v))/{\E\ 
where (...) represents an ensemble average, and 
v={x^ + y^Y'\ 



(2) 



(3) 

With these assumptions all one and two point statistical 
properties can be obtained from W (w) and its deriva- 
tives. 

We take E to be in the far field of a circularly symmet- 
ric planar source S (u), where u measures radial position 
in the source plane. The VanCittert-Zenike theorem [1] 
relates W (v) and S {u) by 

/•OO />CXD 

W {kv) — / uS (u) Jo (kuv) du / uS (u) du, (4) 
Jq Jo 

where n — 2Tr / (AZ), with A the wavelength of light, and 
Z the distance to the remote screen on which E is mea- 
sured. Here and throughout J„ is the Bessel function of 
integer order n. We define 

r = KV, (5) 

and hereafter write W (kv) = W (r). 

S (u) is the intensity distribution (modulus squared) of 
the field leaving the surface of the random medium (sam- 
ple) . We assume this sample to be a random phase plate 
(diffuser) that is deep enough to yield Gaussian statistics, 
but is not so deep as to produce multiple scattering [1]; 
based on this assumption we equate S (u) with the in- 
tensity distribution illuminating the sample, and assume 
that the scattered field retains the state of polarization 
of the incident field. 



A. Simple sources 

We consider initially three different, widely studied 
simple source distributions that individually generate 
normal speckle patterns: a Gaussian (superscript G); a 
disk (superscript D); and a narrow annulus, or ring, (su- 
perscript R); writing: 

(u) = exp [-uV (4p2)] ; (6a) 

^(D) (u) = e (u - p) ; (6b) 

(u) =e6{u-p), (6c) 



throughout Q (w) is the Heaviside step function defined 
by 6 (w < 1) = 1, 6 (w > 1) = 0, and S {w) is the Dirac 
Delta function. 

Although there are no particular problems in generat- 
ing experimentally the simple Gaussian and disk source 
distributions, the ring can only be approximated. This 
can be done by using a narrow, but finite width, annulus. 
Writing S*^^^ (u) for this annulus we have 

S(^) (u) = e (w - p - e/2) - e (u - p + e/2) , (7) 

where p is the mean radius, and e the width, of the an- 
nulus. Upon passing to the limit of zero width while 
noting that 



lim 

£^0 



(u)] /e = 6{u-p) 



(8) 



it becomes evident that for consistency with the dimen- 
sionless Gaussian and disk source functions, the source 



function of the ring must be written as in Eq. ( 6c 



Using Eq. Q the autocorrelation functions W of the 
simple sources in Eqs. (6al-(6cl 



are 



(^pr) 2 Ji (pr) I (pr) . 
(pr) = Jo (pr) , 



(9a) 
(9b) 
(9c) 



whereas the autocorrelation function of the finite width 
annulus, Eq. Q, is [3] 



W(R) (pr) 



: — [(p + e/2) Ji (pr-l-er/2) 
per 

(p-e/2) Ji(pr-er/2)], 



^0 (pr) 



1 



6 



24p 



Ji (pr) + O 




(10a) 



(10b) 



(9c 



As discussed in [3], Eq. (10a I is an excellent approxima- 
tion to Eq. 



for er < 1, so that experiments with 
narrow annuli can be expected to closely match the the- 
oretical results obtained here for a delta function ring. 

The inverse of the parameter p is a characteristic length 
that determines many properties of the speckle field. For 
example, in all cases the number density of extrema, sin- 
gularities, umbilic points, etc. equals cp^/vr, where the 
constant c, which does not differ very markedly from 
unity, depends upon the the form of the autocorrelation 
function and the nature of the critical point. Thus, the 
mean separation between critical points is of order 1 /p, 
the density of level crossings (line length/unit area) is 
proportional to p, and the mean spacing between cross- 
ings along a straight line is proportional 1/p. 
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B. Composite sources 

These important preliminaries completed, we now turn 
to composite source functions that generate speckled 
speckle fields. We build these composite sources out 
of pairwise combinations of simple Gaussian, disk, and 
ring sources. One member of the pair, for which we set 
p — a, has a narrow, but intense, distribution of intensity, 
the other member, for which we set p = b, has a broad 
but weak intensity distribution. 

As a specific example we consider a composite source 
made up out of two Gaussians. We generate this source 
(conceptually) using a single laser beam and an optical 
system that produces two concentric foci at the diffuser 
- an intense tight focus, the a beam, and a weak broad 
focus, the b beam. The resulting setup is shown schemat- 
ically in Fig. |2] Leaving aside the important, practical 
details of the optical system that produces the a and b 
beams, we proceed to analyze the properties of the resul- 
tant speckled speckle fields. 

We start by noting that now there are two different 
characteristic lengths, Aa = 1/a, and A;, = 1/6, with 
Aq ^ Aiy. In Fig. [l] it is the a field that generates 
the large speckle spots, whereas the b field generates the 
small spots that decorate the a field spots. Although 
this result is, obvious, as becomes apparent there are a 
host of nonobvious, often surprising, phenomena that lie 
hidden in these patterns 

Using superscripts Tq and Tf, to denote the source 
type, where Ta,Tb = G, D, or R, we write our composite 
source as 



(11) 



where / is the illuminating intensity (optical power/unit 
area): for a Gaussian / is the peak intensity at the cen- 
ter of the Gaussian; for a disk / is the uniform intensity 
throughout the disk; and for a ring I is the uniform in- 
tensity within the annulus. Inserting Eq. (Ill into Eq. 
Q we obtain for the autocorrelation function of the com- 
posite source 



VF(^°)(ar)-t- (6r) 
1 + x(T„T,) ■ 



(12) 



The dimensionless parameter K^'^'^'^''^ in Eq. (12 1 de- 
pends on the combination of sources that are used, and 
is given in Appendix A. As discussed in this appendix, 
K equals the total optical power in the b source beam di- 
vided by the total optical power in the a source beam, so 
that the composite field autocorrelation function is the 
optical power weighted average of the individual a and b 
autocorrelation functions. 

In addition to W, we will need throughout later sec- 
tions mean square derivatives of the real (imagi- 
nary) parts of the composite optical field, and moments 
of the composite source functions, in both cases 
for n = 1,2,3; r2n is given in Appendix B, M2n in Ap- 
pendix C. We will also need the dimensionless parameter 
r] = rA^lr\ — 1, which is discussed in Appendix D. 



Speckle Field 





Source Plane 



Optical System 



FIG. 2: Apparatus for generating speckled speckle. A laser 
beam enters an optical system that produces two concentric 
beams which illuminate a deep phase plate (the source plane) . 
One beam, labelled a, has a small diameter and a large inten- 
sity, the other, labelled fe, has a large diameter and a small 
intensity. The resulting speckle field is observed in the far 
field of the source plane. Under the paraxial assumption used 
here, the distance Z between the source plane and speckle 
field satisfies Z ^ u,v. The orthogonal xy axis system in 
the speckle field can have an arbitrary orientation, and is not 
shown. The large spots in the speckle field are due to beam a, 
the small spots to beam b. For scalar speckled speckle both 
beams have the same polarization, for vector speckle speckle 
they are orthogonally polarized. Using the above geometry, 
we later illustrate many of the special properties of speckled 
speckle with a computer simulation in which a source that 
consists of 1000 randomly phased points (the Source Plane) 
with a desired composite amplitude distribution (binary com- 
binations of Gaussians, Disks, and Rings) radiates into the far 
field (the Speckle Field). 



III. CRITICAL POINTS OF THE PHASE 



A. Scalar Field Statistics and Morphology 



1. Statistics of vortices and phase extrema 



General results. 

The number density V of vortices [4 — 6], i? of phase 
extrema [7], and S of phase saddles, in normal speckle is 
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extended by inspection to speckled speckle as 

(l 



(13a) 



(TaT(,) 



3/2 



1/2 



(13b) 



where ' is given in Appendix B, M2 and 

m\ " in Appendix C. From the requirement that the 
net Poincare index of all critical points in the plane equals 
+1 follows that the densities of vortices, exrema, and sad- 
dles, are related by 



S=V + E. 



Statistics for single sources. 



(14) 



V and E for a single source, Eqs. (6a I, (6b I, and 



([6c|, are recovered from Eqs. (13a) and (13b I by setting 
K = Q, with results [7]: 
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(15a) 

(15b) 
(15c) 

For a single source we therefore have (£;/]/) = 0.414, 
(£;/y)(D) = 0.0887, {E/V)^^^ = 0, so that E/V is al- 
ways less than unity; for speckled speckle, however, the 
situation is remarkably different. 

Anomalous statistics for compound sources. 

In what follows we always take a < b, and assume that 
p — a/b is small: only in this limit are there two widely 
different length scales in the speckle pattern. We then 
find that for any given small value of p the ratio of phase 
extrema to vortices, {E/V)'-'^'''^'''' = i;(TaT,)/-^(T„T,)^ 
attains a large, maximum value, iE/Vy^^J''\ when 



(E/V) 



(T„Ti,) 
max 



\/3r 



(16) 



with TO2?i given in Eqs. (C5l and (B5) . For example 



for a a Gaussian (G) and b a ring (R) , 



(Elvi^^^ 

V / J max 



(GR) _ ^ 



(G) 



TO 



(R) 



V3r 



,(R) 



12TOf ^TOf) Vp' 



73 

48p' 



,2 ■ 



(17a) 
(17b) 



Comparison of results for single and compound 
sources. 

We compare the above results for speckled speckle 
with those for ordinary speckle. For, say, p = 0.1, 
which should be achievable in practice without undue 



(GG) 



29, 



difliculty, we have: for two Gaussians, {E/V) 
which is 70 times larger than for normal speckle, Eq. 
( |l5a| ; for two disks {E/V)^^^^ « 19, which is 215 larger 
than for normal speckle, Eq. (15b I; and for two rings 
{E/V)^^^^ ~ 14. This latter result is especially striking 
when it is recalled that for normal speckle, Eq. (15c I, 
there are no phase exrema at all. We note that in all 
three of the above cases K^ax = 0.01, which corresponds 
to an intensity ratio h/Ia = 10""* and a total optical 
power ratio Pb/Pa = 0.01. It thus seems likely that the 
major experimental problem in observing the anomalies 
of speckled speckle will be the reduction of parasitic light 
that would tend to swamp the weak b field intensity If,. 

Representative Graphs. 

In Fig. [sjwe display graphs of V, E, and E/V, for two 
Gaussians (GG), two disks (DD), and two rings (RR). 
The graphs for say DR (a = D, 6 = R) or RD (a = R, 6 
= D) can be approximated by dividing the figures along 
the line K = 10~^ and connecting the appropriate com- 
bination of vertically rescaled figure halves; for example, 
the figure for DR is approximated by attaching the left 
half of Fig. ^h) to the right half of Fig. [sjc). 

The graphs in Fig. [3] are representative, and other 
source combinations show similar results: the density of 
vortices {V), Eq. ( |13a[ ), interpolates smoothly from its 
a value (small K) to its b value (large K), as does the 
density (not shown) of saddle points, Eq. (14 1, whereas 
both the density of extrema (E), Eq. ( |13b[ ), and the ratio 
E/V , are anomalous, and reach a maximum for small K. 
The density of these maxima can be made arbitrarily 
large by suitable choice of p and K, Eq. (16 1. 



2. Morphology 



Spatial arrangements of vortices. 

The spatial arrangement of the vortices is also anoma- 
lous. Because a vortex is an absolute zero of the wave- 
function [8] , vortices can appear only where there is com- 
plete destructive interference between the strong a and 
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FIG. 3: Critical point densities, V, E, and E jV , vs. K for 
a = l,ti = 10 (p = 0.1), for (a) two Gaussians (GG), (b) two 
disks (DD), and (c) two rings (RR). Note the logarithmic 
scales. 



weak h fields, so composite field vortices are constrained 
to lie in simultaneously dark regions of the a field, and 
the composite field, speckle patterns. Fig. |4] For very 
small Ihjla the h field is a minor perturbation of the o 
field. In this limit the sparse zeros, and therefore the 
sparse vortices, of the a field maintain their identities 
but are shifted slightly in position. Fig. [4]^ a). For larger 
intensity ratios, however, both a and h field vortices lose 
their identities. 

Reactions between vortices and phase stationary points. 

As lb increases new vortices appear. Conservation 
of charge requires these vortices to nucleate as posi- 
tive/negative pairs. Pair production by itself, however, 
violates conservation of the Poincare index, here Jp, be- 
cause Ip — +1 for vortices independent of their charge. 
Vortex pair production must therefore be accompanied 
by another process, and in general there are two possi- 
bilities: 



I. A vortex pair and a pair of saddle points (each one 
of which has Ip = —1) are generated [9]. 

II. Two extrema (one minimum and one maximum, 
each of which has Ip = +1) collide and annihilate, and 
a vortex pair emerges from the collision [10]. 

For very small Ij, there are, in general, fewer extrema 
than vortices, so in this regime mechanism II cannot be 
operative (recall that for a ring there are no phase ex- 
trema at all), leaving mechanism I as the only alterna- 
tive. In Fig. [3] this regime corresponds to the region of 
say K < 10~^. But as can be seen from this figure, for 




FIG. 4: Speckled speckle phase vortices for Gaussian a and b 
fields (GG) with p = a/b = 0.1. Shown is a typical region 
of the intensity speckle pattern with vortices superimposed, 
where filled white circles with black rims (filled black circles 
with white rims) represent positive (negative) vortices, (a) 
h/Ia = 2 X 10"^ {K = 0.002). The large circles show the po- 
sitions of the vortices of the strong a field (It, — 0) , the small 
circles show the vortices of the composite field. As expected, 
for the small lb/ la ratio used here, the composite field vortices 
cluster near the positions of the sparse a field vortices, which 
are the darkest regions of the speckle pattern. As h/Ia is 
increased from zero, and as is illustrated here, initially the a 
field vortices undergo small displacements but maintain their 
identities, while a small number of new positive/negative vor- 
tex pairs nucleate (here two pairs, arrows), (b) The field 
in (a) for lb/ la = lO^"' {K — 0.1). As lb/ la increases the 
immber of composite field vortices increases rapidly (Fig. [3|, 
and in addition to vortex nucleation vortex annihilation also 
occurs, leading to a loss of identity of the a field vortices. For 
the moderately small h/Ia ratio used here the composite field 
vortices no longer need cluster around the positions of the a 
field vortices; they are, however, still constrained, to lie in the 
dark regions surrounding the large o field speckle spots. For 
still larger K vortices begin to fill the field more uniformly, 
and in the limit of large K the dense vortex structure is that 
of the 6 field. 
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FIG. 5: Phase maps showing vortex nucleation for smaU and 
intermediate K for Gaussian a and b source fields (GG) with 
a/6 = 0.1. The phase increases to 27r black to white. Filled 
black (white) circles represent negative (positive) vortices, 
whereas black squares mark saddle points. Thin curves are 
ordinary contours, thick curves are the special contours (bi- 
furcation lines) that pass through the saddle points. (a),(b) 
Vortex nucleation for very small K. (a) K = 0.00164. (b) 
K = 0.0017. In this small K region mechanism I of the 
text is dominant, and vortex pairs nucleate together with a 
pair of saddle points. The resulting box like structure, first 
described in [9], is the one usually observed following vortex 
pair production in normal speckle, (c)-(f) Intermediate K 
regime. In (c), (d), (e), and (f), K = 0.024,0.0275,0.0291, 
and 0.0325, respectively. In this region mechanism II of the 
text dominates. Specifically: in (c) there are no vortices, only 
phase extrema (two maxima and one minimum); in (d) the 
upper minimum and maximum approach one another, while 
in an unrelated process a new saddle point and a new min- 
imum nucleate; in (e) the upper pair of extrema collide and 
generate a positive/negative vortex pair; and in (f) the lower 
pair of extrema collide and generate a second vortex pair. 



larger K the number of extrema exceeds the number of 
vortices, thereby permitting mechanism II to become op- 
erative. These two reaction mechanisms are illustrated 
in Fig. |5]for the representative case of GG for small and 
intermediate K, and in Fig. [6]for RR for larger K; as can 
be seen from these figures, for sufficiently small (large) 
K mechanism I (mechanism II) dominates. 



B. Vector Field C points and Extrema 

1. Introductory remarks 

The state of polarization of generic vector fields is el- 
liptical. In the paraxial fields of interest here the (always 
planar [11]) polarization ellipses lie in planes, here the 
xy-plane, oriented normal to the direction of propaga- 
tion, the Z-axis. The point singularities of such fields 



in the plane (in space) are C points (form continuous C 
lines), which are points (lines) of circular polarization em- 
bedded in a field of elliptical polarization [12]. Because 
the azimuthal orientation of a circle (the C circle) is un- 
defined, the surrounding ellipses rotate about the point, 
generically with rotation (winding) angle ±7r, and topo- 
logical charge (winding number) ±i for positive/negative 
points [12]. 

Necessarily associated with C points in the planes of 
generic paraxial fields are azimuthal stationary points - 
extrema (maxima and minima]), and saddle points [13]. 
Unlike C points which are easily identified visually as the 
centers around which the surrounding ellipses rotate, az- 
imuthal stationary points are subtle features that may 
be difficult to discern in traditional maps that show po- 
larization ellipses. A highly useful representation that 
shows both types of critical points clearly is the phase 
(argument) <i>i2 of the complex Stokes field [14] 

Si2 = Si + iS2 (18) 

In this representation C points are vortices with charges 
±1, and azimuthal stationary points are phase maxima, 
minima, and saddle points. In what follows we use this 
representation exclusively. 

Charge conservation requires that C points nucleate or 
annihilate as positive/negative pairs, whereas conserva- 
tion of the Poincare index Ip, which is -t-1 for C points 
of either sign, requires that pair production be accompa- 




FIG. 6: Phase maps showing vortex nucleation for larger 
values of K for both a and b source fields rings (RR) with 
a/b = 0.1. Symbols are as in Fig. [5] (a) = 0.5. (b) 
K = 2. In (a), which corresponds approximately to the re- 
gion of the fiat maximum in the number of extrema E in Fig. 
|3j;, there are three vortex pairs and numerous maxima and 
minima. In (b), which corresponds to the onset of rapid de- 
cline in the number of extrema in Fig. [Sj;, all ma:xima and 
minima in (a) have collided, annihilated, and been replaced 
by vortex pairs. 
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nied by either the production of two azimuthal saddles, 
for which Ip = —1, or the disappearance of two extrema, 
for which Ip = +1. 

Like ordinary circularly polarized light, C points are 
either right handed or left handed, depending upon the 
direction of rotation (clockwise or counterclockwise) of 
the electric vector of the light as it traces out the C circle. 
By continuity, the polarization ellipses in the immediate 
vicinity of a C point have the handedness of the point, 
so that overall the ellipse field is divided into right and 
left handed regions. Near the boundary between regions 
of opposite handedness the polarization ellipses narrow, 
and on the boundary itself the eUipses collapse into hues 
and the polarization is linear; these lines are called L lines 
[12]. 

All of the above holds for both ordinary, and speckled 
speckle, vector fields. 

We now turn to a discussion of the number densities of 
C points and their associated extrema in vector speckled 
speckle. Here the a and b source fields have orthogonal 
polarizations, and we consider two different cases: in the 
first the a (b) source field is linearly polarized along the 
X-axis (y-axis); in the second the a (b) source field has 
right (left) circular polarization. As becomes apparent, 
these two different configurations produce very different 
results. 



Linearly polarized source fields: statistics and ellipse 
field morphology 



a. C point densities. 

In what follows we assume, as already stated, that 
speckle field {Ey) maintains the state of linear po- 
larization of its a (b) source. We begin by noting that at 
a C point the ampHtudes of and Ey must be equal, 
so like phase vortices, C points of the composite field are 
found in the dark regions of the a field speckle pattern 
[2]. But unlike vortices, these regions need not also be 
dark regions of the composite field speckle pattern, be- 
cause unlike a vortex, a C point is not a zero of the vector 
field. 

In considering number densities it is convenient to de- 
compose the speckle field. Fig. [2] into right (R) and left 
(L) circularly polarized components, Ep and E^. The 
reason is that Ep {El) vanishes at a right handed (left 
handed) C point, and since the single polarization com- 
ponent fields Ep and El are scalar fields, Eq. (13a I is 
applicable to each component. From this follows that in 
normal speckle fields the number density of C points is 
twice the density of phase vortices [12]. 

Eq. ( 13a I is also applicable to C points in vector speck- 
led speckle, but to use this result we need derivatives 

^(ToT,) 



for the circularly polarized components. Using 
the results in Appendix B, we can obtain these deriva- 
fives from the autocorrelation functions VF]j £ of Ep 



and El', these autocorrelation functions are calculated 
below. 

We start by writing the circularly polarized compo- 
nents Ep and El in terms of the Cartesian components 
Ex and Ey of the speckle field as 



(19) 



Ep,L = (E, T ^Ey) /V2, 



where E^ (Ey) arises from the x-axis (y-axis) linearly 
polarized a (6) source field. Temporarily suppressing for 
notational convenience superscripts TaTf,, Wp^L is 



Wp,L (r) = {E*p^L (0) Ep^L (r)) / {Ip^i 



(20) 



where the average intensity {Ip l) = 
{E*^,^{Q)Ep^l{Q)) = \{{h) + {Iy)), with (/,) and 
{ly) the average intensities of the Cartesian compo- 
nents. Using Eq. ([l9) together with {El (0) Ey {r)) = 

{e; (0) Ex (r)) - (irfo) Ey (0)) = {e; (O) E^ (O)) = o, 

we have 



R 



) (r) = VFf "^'^ (r) = M^(T„T,) ^ (21) 



where l^C^"'^'') (r) is given in Eq. (12 1, and K is given 
in Eq. (A3 1. In obtaining this result we use the 
fact that {Ix} {{ly)) is proportional to the total opti- 
cal power in the source field a {b ) component, so that 
{ly) I {Ix) ^ K. Thus, Eq. ([I3a| and curves V in Fig. 
[3] multiplied by a factor of two, also describe C points in 
vector speckled speckle. 

b. L line densities. 

The density (line length/unit area) of L lines is pro- 
portional to the square root of the density of C points 
[15]. Because the density of C points, 2V , interpolates 
smoothly between its small K {a field) and large K {b 
field) limits, so does the L line density, and L lines in 
speckled speckle do not show anomalies. 

c. Densities of azimuthal extrema. 
K and p = a/b both small. 



Although Eq . (|13a l is applicable to speckled speckle 
C points, Eq. ( 13b I is not necessarily applicable to the 



corresponding azimuthal stationary points; the reason is 
that unlike the case of C points, there is no a priori sim- 
ple relationship between vector field azimuthal station- 
ary points and the stationary points of the corresponding 
scalar field components. At present there is no exact the- 
ory for the azimuthal densities; we are, however, able to 
estimate these densities when both K and a/b are small. 
Returning to Eq. ( 19 1 we write 



^2Ep^L = (1 T iEy/Ex) 



(22) 



where the motivation for this form is the fact that when 
K is small so also is Ey/E^. Writing the various field 
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components in terms of amplitudes A and phases cp or ^, 



Ex = exp {iipx) , 
Ey = Ay exp (iipy) , 
Er = Ar exp {i<^R) , 
El = ALexp{i^L) , 



(23a) 
(23b) 
(23c) 
(23d) 



and noting that Stokes parameters 5*1 and ^2 can be 
written 



51 = 2Re{E*REL), 

52 = 2Im iE*REL), 



we have from Eq. ( 18 1 



Si2 — 2E*j^El, 



from which follows 



$1 



(24a) 
(24b) 



(25) 



(26) 



Returning to Eq. ( 22 1 we have to third order in the 
small quantity Ey/E^, 

1 ± tEy/E, « exp {{Ey/Exf/2) 



Writing 



X exp Ti Ey/E, - (Ey/E^y /3 . (27) 



Ey/E^ = {Ay/ Ax) [cos {iPy - If x) + 1 COS {iPy ^ (fix)] 



we have 

^R,L ^ ^PxT {Ay /Ax) cos {(py - ipx) 
+ \{Ay/Ax)'^sin (2 {ipy - tpx)) 

±\{Ay/Ax)^COS^{iPy-iPx) 

=F {Ay/Ax)^ sin {ipy - ipx) sin (2 ((^^^ 
so that to leading order in Ay /Ax, 

'^12 



(28) 



^x))- (29) 



~2 {Ay /Ax) cos {ipy - (fix) (30a) 
2 ($L -</=.) (30b) 
-2{^R-iPx). (30c) 



Eqs. (30 1 make a number of possibly unexpected predic- 



tions that are listed below. 

Density of azimuthal extrema for K equal to its value 
at the peak of curve E in Fig. [5[ 

Far from its sparse vortices the phase ipx of the slowly 
varying a field changes slowly, so that in regions where 
'^R,L vary rapidly, tpx in Eqs. (30b I and (30c I can be 
treated locally as an unimportant, nearly flat background 
that can be eliminated by a local redefinition of the zero 



of phase. For sufficiently small p = a/h these flat regions 
occupy most of the wavefleld area and we can write 



$1 



2$, 



(31) 



In the region of the peak of curve E in Fig. [3] where 
there is a high density of extrema, $_r,l vary rapidly. We 
therefore conclude that near this peak the density of az- 
imuthal extrema very nearly equals the density of phase 
extrema E, so that just like for scalar speckle speckle, 
E /V for vector speckle speckle can attain very large val- 
ues. 

In Fig. [7] we compare (px, '^12, '^l, and —'^r for the 
case of two disks (DD) for the small, but still experimen- 
tally possible, ratio a/h = 0.01, for K = 1 >(. 10~^, which 
corresponds to the maximum of E /V for the above a/h 
ratio (here {E/V)^^^ = 1924). These parameters are fa- 
vorable ones, and for larger a/h ratios, say a/b = 0.1, the 
very close agreement with the predictions of Eqs. (31 1 



degrades, because now Lpx can no longer be considered to 
be locally flat over most of the wavefleld area; the over- 
all conclusion, however, remains, and in general E/V is 
anomalously large for vector speckle speckle. 

Density of azimuthal extrema for K less than its value 
at the peak of curve E in Fig. [?[ 

For small ratios a/b for all values of K less than Kq, 
the value of K at the peak of the E curve, the stationary 
point structure of <i>i2 remains very nearly independent of 
K . This follows from the fact that not only the phase ipx 
of the slowly varying a fleld hardly changes, but that also 
Ax, the amplitude of the fleld, changes slowly, and can 
therefore locally be taken to be constant. Eliminating 



FIG. 7: Azimuthal extrema in vector speckled speckle for 
a/h — 0.01, and K — 10"''. The amplitude in these contour 
maps increases black to white, so that high maxima (deep 
minima) are colored white (black), (a) ipx- The variation 
in ipx is less than 7% of the variation of "l>i2 over this same 
region, so to first order ipx can be taken to be constant in Eqs. 
(301. (b) $12. (c) — $_R. The minus sign turns maxima into 
minima and vice versa, facilitating comparison with "l>i2 and 
$L. (d) ^L- Note the close similarity of the three fields $12, 
$L, and $L, in full accord with Eqs. (31 1. 
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ipx by redefining the local zero of phase, and replacing 
Ax by a local constant c, we have 



'CAy COS {iPy) . 



(32) 
and the 



Now, for a given value of p, the structure of Ay, 
structure and magnitude of ipy , are independent of la and 
/{,, and are therefore independent of K; Eq. (32 1 there- 
fore implies that although the magnitude of <f>i2 changes 
with K, its structure remains invariant. This is illus- 
trated qualitatively in Fig. [sjja). 

In Fig. ^h) we quantify the relationship between $12 
and ^L- Shown is the cross correlation coefficient C 
between $12 (AT) and ^l{Ko), where C = C {u,v) for 
two real variables u, v is 



C{u,v) 



{iu-{u)){v-{v))) 



(33) 



Here (...) represents spatial averages over the phase fields. 
As can be seen, $12 and are nearly identical (C — 
0.983) up to and including the point at which E passes 
through its maximum, showing that in vector speckled 
speckle the density of azimuthal extrema all throughout 
the small K region substantially equals the maximum 
value of E for scalar speckled speckle. 

An interesting, possibly counterintuitive, feature of 
Fig. ^h) is the significant degree of correlation (C — 0.8) 
between <i>i2 and <&l in the large K region where most 
extrema have been converted to C points. This corre- 
lation reflects the fact that locally the phase field before 
and after collisions of extrema maintains more-or-less the 
same value. This behavior is most clearly seen in Fig. [6] 
where dark (light) regions in Fig. [6j^a) for the most part 
remain dark (light) in Fig. [6]^b). 

Generation of C points as K increases from zero. 

As noted above, the density of C points increases with 
K in accord with curves V (multiplied by two) in Fig. [S] 
For scalar speckled speckle vortex generation in the small 
K region in which extrema are sparse occurs via mecha- 
nism I in which a new vortex pair appears together with 
a pair of saddle points. Figs. [5][a),(b), whereas in the 
region of larger K in which extrema are plentiful vortex 
pairs are generated via mechanism II in which the colli- 
sion and annihilation of a pair of extrema, one maximum 
and one minimum, generates the vortices. Figs. |5|^c)-(f) 
and Fig. |6] For C points, however, extrema are abun- 
dant for both small and large K , and it is mechanism 
II, i.e. collision of azimuthal extrema, that dominates C 
point production for all K. 

d. Decrease in the number of azimuthal extrema for large 
K. 

In the limit of large K where E^/Ey is small, the ana- 
log of Eq. (l22| is 



and the analogs of Eqs. (|30| are 

$12 w -2 {Ax/ Ay) cos {(fix 
«2(ci>i-^y), 



(35a) 
(35b) 
(35c) 



For large K, in scalar fields <i>i, <^r, and ipy, E decreases 
from its large maximum to its small single component 
value (Fig. [3|. Eqs. (35) therefore imply that also 



the large density of extrema present in vector speckled 
speckle for small K decreases for large K. For scalar 
fields the reduction in the number of extrema results from 
collisions of maxima and minima that produce vortex 
pairs. Figs. [5] and [6] as is illustrated in Fig. [sj^c), the 
same is true for vector speckled speckle. 

e. C points for vanishingly small values of either the a 
or the b field. 

As K increases the phase structure stabilizes asymp- 
totically at a mix of C points and extrema, with C points 
outnumbering extrema by perhaps a factor of 3 — 5. But 



if for, say, \Ex 



10 



-25 I 



Ey\ {K ^ 10^4), how, for all in- 



tents and purposes, can the wavefield be anything other 
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V2Er^l = T^Ey{l±iEx/Ey), 



FIG. 8; Azimuthal extrema for two disks (DD) for a/b — 0.01, 
as in Fig. [t] (a) $12 for K = lO"'^". Note that the station- 
ary point structure here is very nearly the same as in Fig. 
[rj^b) where K = 10~* is fifty orders of magnitude larger, (b) 
Correlation coefficient C = C ($l {Kq) , $12 (K)), Eq. 
"I>L is the phase of the left circularly polarized vector field 
component, "l>i2 is the Stokes phase of the combined vector 
field. Ko — 10~^ is the value of K at which the scaled number 
density E/Emax of $_l extrema, dotted curve, attains its max- 
imum value. Curve E/Emax is plotted as l + logj,, (-E/iSmax), 
so that where C = 0.98 (C = 0.8), E/E^^^ = 0.95 (£/£max = 
0.63). Note the extreme range of K, 10"^^ <K < 10^". (c) 
$12 for K = 1. As can be seen, when K starts to become 
large the azimuthal maxima and minima in (a) start to col- 
lide and annihilate, generating positive/negative C point pairs 
(filled white/black circles). 



(34) 
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than linearly polarized? The answer [12, sect. 12.5] is 
that Ey goes continuously to zero at its vortices, so that 
surrounding each Ey vortex, which is an ly — \Ey\^ zero 
minimum, there is for every K an elliptical contour on 
which \Ea: \ = \Ey\. The size of this contour shrinks with 
increasing K, so for sufficiently large K the phase (px of 
Ex on the contour can be taken to be constant. Because 
the phase field surrounding the Ey vortex winds through 
27r there are always two points on the contour (sepa- 
rated by 180°) where \(py — (px\ = 7r/2; at these points 
a positive/negative pair of C points appears. Similar 
considerations apply to the limit of very small K, where 
the C points are located near the vortices of E^ ■ 

f. Comparison of results for scalar and vector fields. 

We conclude this subsection with a summary of the 
similarities and differences in the number densities of vor- 
tices, C points, and extrema, in scalar and in linearly po- 
larized speckled speckle: (i) For both cases the density 
of singularities - vortices for the scalar case, C points for 
the linearly polarized vector case - interpolates smoothly 
with K from its small a field value to its large b field 
value: Fig. |3] curve V for vortices; this same curve mul- 
tiplied by two for C points, {ii) For the scalar case the 
density of extrema, and the ratio of extrema to vortices, 
increase from their a field values to large maxima at small 
K, and then decrease for large K to their b field values: 
Fig. [sj Eqs. (16 1, {in) For the vector case the density 
of extrema, and the ratio of extrema to C points, are 
large for all small K, and decrease for large K. For suf- 
ficiently small a/b the large azimuthal extrema density 
for small K substantially equals the maximum density of 
phase extrema for the scalar case, whereas the large ratio 
of azimuthal extrema to C points equals half the maxi- 
mum ratio of phase extrema to vortices for the scalar 
case, (iv) In the scalar case generation of vortex pairs 
is accompanied by nucleation (annihilation) of a pair of 
saddle points (extrema) for small (large) K. For the 
vector case C point generation in speckled speckle occurs 
through the annihilation of paired azimuthal extrema for 
all K. (v) The reduction in the number of extrema for 
large K results from collisions of maxima and minima to 
produce vortices in the scalar case, C points in the vector 
case. 



which implies the following: (i) The number densities, 
and the spatial distributions of all critical points - C 
points, azimuthal extrema, and azimuthal saddle points - 
are independent of K, which depends on the amplitudes 
of the fields, but not on their phases. (ii) For small 
values of a/b for which the slowly varying a field phase 
ipa can, over most of the wavefield area, be treated lo- 
cally as a constant, the critical point structure of $12 
is essentially the same as that of (pb- In this limit the 
critical point number densities are those of scalar field 6, 
and are given by Eqs. ( 15 1 with p — b. These conclu- 



sions are fully verified by our computer simulations (not 
shown). Thus, in contrast to linearly polarized speckled 
speckle with its rich set of anomalies, circularly polarized 
speckled speckle appears to offer little that is of special 
interest. 



C. Zero Level Crossings of Phase Derivatives 



Introductory remarks. 

Nodal domains and the zero level curves that define 
them are of considerable current interest because they 
determine many important properties of random fields 
[1]. At present there is little information available on 
the nodal domains of ordinary vector speckle; this pre- 
vents us from using our present methods to calculate 
these domains for vector speckled speckle, and below we 
consider nodal domains of the phase only of scalar speck- 
led speckle. 

Densities of zero crossings and densities of intersec- 
tions of these crossings with a straight line. 

Phase extrema, saddles, and vortices [16], are all lo- 
cated at intersections of the zero level crossings and 
Zy of the first order phase derivatives (fix and tfy. It is 
therefore reasonable to expect that like the number den- 
sity of extrema, E, which shows an anomalous peak in 
Fig. [3]at iiT ~ 0.13, also the line densities Dx — Dy = D 
of Zj. and Zy would show similar anomalous peaks. 

For many types of random curves in an isotropic field 
there is a simple relationship [15] between 1) and the 
density N of intersections of these curves with a straight 
line, say the x-axis. 



3. Circularly polarized source fields 

We take the a (6) field to be right (left) circularly po- 
larized, and write the circular components in the speckle 
plane as 



Er = Aa exp{iipa) , 
El = Ab exp {iipb) . 



Eq. ( 26 1 then becomes 

^'12 



(36a) 
(36b) 



(37) 



V ■ 



-N. 



(38) 



This relationship, however, does not hold for level curves 
of derivatives because, as discussed below for the phase, 
differentiation destroys the required isotropy. 



1. Densities of intersections of zero crossings of phase 
derivatives with the x-axis. 



General results. 
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Previously, we calculated for normal speckle the quan- 
tities Nx{(px,L) {Nx {(fiy, L)), which are the density of 
intersections of crossings at level L of phase derivatives 
fx (Vj/) with the a;-axis [17]. Extending by inspection 
these results to the zero crossings Zx and Zy of speckled 
speckle, we have 



This reversal is a general feature, and 
occurs also for two Gaussians (GG) and two rings (RR), 
both of which yield graphs that are quite similar to Fig. 

Explanation of anomalous reversals. 



.(T^Tt) 



.(T„T6 



1/2 



7T 

where rj"^""^""^ is given in Appendix B, TyC^""'"'') in Ap- 
pendix D. 

Now, in the isotropic random field of interest here, 
on average the length/unit area of Z^, and Zy must be 
the same, therefore, if Eq. (38 1 were to hold we would 
have N^^"'^"'' {(px,0) = N^^"^'' {ipy,0); but this is not 
the case, and instead [17] 



xE 




.(ToT,) 



<2, 



xE 



l(TaTi) 



V 



(40) 



where E is the complete elliptic integral of the second 
kind. 

Results for a single source. 



Eqs. (39) and ( |40[ ) with single superscripts hold, of 
course, for ordinary speckle. For a Gaussian (G), disk 
(D), and ring (R), 77^^) = 2, 77^^) = 1, t]'-^'^ = 1/2, and 



i.e. < N, 



(T) 



N!i':J/Ni%^ = 0.900, 
Ni^^/Ni'^y' = 0.699, 
Ni'l^/NiT ^ 0-524, 
in all three cases. 



(41a) 
(41b) 
(41c) 



Anomalous reversals for compound sources. 



From Eqs. (|39|, (|40|), and (|B4|, it follows that b'^Nx^x 
and b^^Nx^y are functions only of p — a/b and K. Ac- 
cordingly, in Fig. 9|a) we plot b^^N^^J^^ and b^^Nx^J^'^ 
vs. K for two disks (DD) for p = 0.1. As can be seen, 
for very small and very large K (the regions labelled 'b' 

i n th e figure) Nx 



DD) 



(DD) 



in agreement with Eq. 
(41b I. In the intermediate region lO""* < K < 1 (region 



'c'), however, this situation reverses, and anomalously 



The anomalous Nx^x/Nx^y ratio, as well as the normal 
ratio, find their explanation in the structure of the phase 
field. For normal speckle, and in regions 'b' of Fig. [9](a) 
for speckled speckle, vortices outnumber extrema. As 
discussed in [16], in the immediate vicinity of a vortex 
Zx (Zy) is necessarily strictly parallel to the a;-axis (y- 
axis). Fig. [9]jb). As a result the number of Zx (Zy) that 
cross the x-axis (y-axis) is reduced (increased), resulting 
in Nx.x < Nx.y. But in region 'c', extrema outnumber 
vortices, and because, as discussed in [18], at extrema 
there is a strong tendency for Zx {Zy) to be biased per- 
pendicular (parallel) to the x-axis. Fig. ^c), the number 
of Zx (Zy) that cross the x-axis is increased (reduced), 
resulting in Nx,x > Nx,y. 
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FIG. 9: (a) Densities of zero level crossing with the a: -axi s 
for two disks (DD). Thin curve, b~^jVi°°' i'fix,0), Eq. ([39 1 
thick curve, b'^Nf°°^ {v>v,0), Eq 
(b),(c) Zero crossings Zx of ip, 



(401 



( z 



Here p — a/b — 0.1. 
of (py), shown as thick 



white (black) lines, superimposed on a map of (b) a phase 
vortex, (c) a phase maximum. In both (b) and (c) the phase 
increases — 2tt black to white. For vortices, Zx (Zy) is 
always strictly parallel (perpendicular) to the a;-axis. For 
extrema, the orientation of these zero crossings depends on 
the anisotropy and orientation of the elliptical contours in the 
immediate vicinity of the extremum. For circular contours, 
and for elliptical contours whose principal axes parallel the 
xy-axis system, Zx (Zy) is parallel to the y-axis (a;-axis). For 
other cases there is a strong bias towards an orientation in 
which the absolute value of the angle that Zx (Zy) makes 
with the a;-axis is greater than (less than) 45°. 
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Relationship between densities of phase critical points 
and zero crossings of phase derivatives. 



Eq. (38 1 does not strictly apply to Tj^ and Zj, sepa- 
rately, nonetheless, because the differences in Nx.x and 
Nx^y are not severe, and in light of the explanation given 
above for this difference, we can reasonably assume that 
Eq. ( 38 1 does apply, at least approximately, to the aver- 



age X-axis crossing density 



T„T6 



(42) 

Because all phase critical points - extrema, saddles, 
and vortices - lie on intersections of Z^, and Zj,, one 
can reasonably expect some relationship to exist be- 
tween the densities of these zero lines and the density 
of critical points. We therefore attempt to find a re- 
lationship between the t otal phase cri tical point density 

E + S = 2S 



V + E + S = 2S,¥.qs 
Noting that V, S, and E 



(13) and (14 1, and 
Tor normal s 



m 

peckle are propor- 



tional to p^, with p the characteristic inverse length scale 
of the autocorrelation function, whereas (N^) is propor- 
tional to p, we make the simplest possible ansatz 



,JiV + E + S) 



(T,Tb 



,(T„ 



1 

2 

5-1 /Ar(T»T. 



(43) 



where c^"^°' and c'-"'"''^ are constants calculated from the 



single component analog of Eq. ( 43 ) (without the factor 
of |) for normal a and b speckle fields using the known 
values of V, E, S, Nx,x, and Nx,y for these fields. In Fig. 
10 we plot the R.H.S. of Eq. (43) as continuous curves 



on which are superimposed small circles representing the 
L.H.S of this equation. Shown are representative results 
for the case of two Gaussians (GG) and the case where a 
is a ring and 5 is a disk (RD). Other source combinations 
that are not shown yield similar results. Considering the 
wide range of densities (10^^ - 1), the broad parameter 
range for K (IQ-^^ - 10^), for p (0.001 - 0.1), and for 
77 (0.5 - 10^, Appendix D), and the fact that there are 
no adjustable parameters, the agreement with Eq. (43) 
seems satisfactory. 



IV. CRITICAL POINTS AND ZERO LEVEL 
CROSSINGS OF THE REAL AND IMAGINARY 
PARTS OF THE WAVEFUNCTION 

The real (imaginary) part of the wavefunction of a 
speckle pattern represents an archetypal Gaussian ran- 
dom surface. The critical points that define the major 
features of such a surface are its extrema (maxima plus 
minima), saddle points, and umbilic points, the latter 
being points of isotropic curvature; these are frequently 
located on the sloping sides of extrema. Here we study 
the number densities of these features, and the densities 



10 



10 



10 



E l 1 1 ll_,J^-,-„-JU>-rA 


- (a) ^-CP-^ 


jf GG 


- 

-= 










p = 0.001 

p = 0.01 








(¥^111 


1 1 1 





10 r 



V 



(b) 



10' 



10' 




p = 


0.001 




0.01 




0.1 



-12 -10 -8 -6 -4 -2 2 4 

10 10 10 10 10 10 10 10 10 
K 



FIG. 10: Test of Eq. (431. The average zero level crossing 
densities, {^x), are shown by continuous curves, whereas 
the square roots of the total density of all phase critical points, 
\/V + S + E, are shown by small white circles. In (a) the a 
and b source fields are both Gaussians (GG). In (b) the a 
source field is a ring and the b source field is a disk (RD). 



of the associated level crossings, for speckled speckle. In 
what follows we consider the real {TZ) part of the wave- 
function, but all results hold equally well for the imagi- 
nary (X) part. 



Critical Points 



General results. 



The number density £ of extrema [19], and the den- 
sity U of umbihcs [20] , in normal speckle is extended by 
inspection to speckled speckle as 



1 rf"^'^ 



1 M 



(44) 



(45) 



Results for single sources. 
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For normal speckle [19,20] Eqs. and (|45| with 

= 0, yield 



— ^(D) 

— - 



37r 

= 18^^' 
1 , 



P 



(46a) 

(46b) 
(46c) 



where 1 /p is the characteristic length scale in the speckle 
pattern, Eqs. The number density ratios U/£ for 

a Gaussian (G), disk (D), and ring (R), are 2.6, 1.9, 
and 1.7, respectively. Thus, in normal speckle each ex- 
tremum (maximum or minimum) is decorated by approx- 
imately one umbilic point; for speckled speckle, however, 
the number of umbilic points per extremum can be arbi- 
trarily large, as we now show. 

Anomalous ratio of umbilic points to extrema in speck- 
led speckle. 



Using Eqs. 



(|44|) and (|45| we find that for any 

max 1 



given small value of p — a/b, when K — JC 

= {U^/E™)^^^ attains a large 
maximum value, where to leading order, 



^4 



{U/£) 



max 




(47a) 
(47b) 



As an exam- 

16" ' 



and 



with C2„, n = 1,2,3, given in Eqs. (|B5). 
pie, for a a ring and h a Gaussian, 
{U/£)^^2^ = 3^/3/p2. 

For T = G, D, R, [U / / {U / E)'^''\= l/ (4p2). 
For p = 0.1, for example, the ratio of umbilic points to 
extrema is therefore enhanced by a factor of 25. 

In Fig. 11 'a) we plot U, £, and U/£ for two Gaussians 
(GG) for p = 0.1. As can be seen, both U and £ ex- 
trapolate smoothly between their limiting a and b field 
values, whereas U/£ shows the expected large peak at 
K = 10~^. These curves are representative, and similar 
results are obtained for other GDR combinations. In 
Figs. [Tl|b) and (c) we show the real part TZ of the wave- 
field with its umbilic points superimposed. These points 
were located using the fact that they lie at the intersec- 
tions of the zero [20] lines of the Stokes-like parameters 
El and T.2, where 



5ji — T^xx 

S2 — 27^2;ij 



n 



(48a) 
(48b) 



and the subscripts imply differentiation with respect to 
X and y. 
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FIG. 11: Extrema and umbilic points for two Gaussians (GG) 
for p = ajh — 0.1. (a) Densities U of umbilic points, £ of ex- 
trema, and their ratio U /£. (b),(c) Gray scale coded contour 
maps of the real part of the wavefunction TZ with superim- 
posed umbilic points (small black circles). TZ increases black 
to white; the small white diamond marks a saddle point, (b) 
h = {K = 0). (c) h/Ia = 10"^ {K = 10-^). Note 
that the underlying stationary point structure (extrema and 
saddle point) in (b) and (c) are very nearly the same. This is 
typical, and even for K ^ 1 the underlying stationary point 
structure of the a field remains visible The 152 umbilic points 
in (c), which corresponds to the peak oiU/£ in (a), are con- 
sistent with the number expected for a region that contains a 
central extremum (here a maximum), and parts of the sloping 
sides of two other extrema (here minima located in the upper 
left and lower right corners of the figure). The hyperbolic 
curvature of the saddle point tends to inhibit formation of 
umbilic points because the principal radii of curvature of the 
saddle have opposite signs, whereas at an umbilic point both 
radii must be equal. 



B. Zero Level Crossings 

General results. 

Level crossing densities for TZ, and its first, TZx, TZy, and 
second, TZxx^ TZxyjIZyy, derivatives have long been known 
for normal speckle [17,21], and we obtain by inspection 
for level zero for speckled speckle. 

Ml,™ ^ AAjT^T.) 0) ^ -Jr^™, (49a) 

TT ' 

(49b) 
(49c) 
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and 

j^!,™ (7^„o) = ^^fi™ (7^.,o) , (soa) 

Af^™ (7^.„0) = yfA^iT-T.) (7^^^,o) , (50b) 

J^i™ (7^,„0) = /^AAi^^^^) (7^..,0) . (50c) 
For the ratios of zero crossing densities we have 

_^(T.T.) 



*r(T„T.) 

■'^1,0 



_^(T„T.) 

where 7y('^<»'^!') is given in Appendix D, and 

■'^2,1 



(T„Ti 



Ml 



Results for single sources. 



33/4 V £:(T„Ti,) 



(51) 



(52) 



In normal speckle fields zero crossings of R, Rx, and 
Rxx have similar densities. Setting K = (Ii, = 0) in 
Eqs. (|51|) and (52 1, we have in obvious notation 



■^1.0^ = V3 « 1.73, Af^f = v/573 ~ 1-29, (53a) 
AA^fD) = \/2 w 1.41, Af^f = 75/2 w 1.12, (53b) 
7V[o^ = V372 « 1-22, M^^^^ = VTo/3 w 1.05, (53c) 



For speckled speckle, however, these ratios can be arbi- 
trarily large, as we now show. 

Anomalous ratios for compound sources. 



A/^Q-^'^ Eq. (511, and Af^^^^^'K Eq. (52 1, can be 



made arbitrarily large by a suitable choice of ^TTaTb) g^^^ 
p. For Ml Q we find that for any given small value 
of p, when A'(TaT5) ^ (/cj^o"^"') , M^l'^'^ attains a 

V ■ / max ' 

-^10'' ^ ) ' where to leading order 

' / max 

in p, 



^ 



1 

2^\ 



„(Ta) (Tb)- 



(54a) 
(54b) 



Similarly, for A/'2^°'^'''' we find 



^4 

L4 



1 

2^\ 



„(T.) (Tb) 

--2 ""e 



(55a) 
(55b) 



As representative examples we plot in Fig. 12 M^ q 

and A/'2^*^'' for p = 0.01. As can be seen, both densities 
show large enhancements at their peaks over their normal 
speckle values. Below we inquire as to the sources of 
these enhancements. 



(TaTt) 



The zero crossings of TZ, which obey Eq. (38), divide 
the field into nodal domains of opposite sign. In normal 
speckle maxima (minima) mostly occupy positive (neg- 
ative) domains. Mx (TZ, 0) therefore roughly scales with 
the square root of the density of extrema. The same 
is true for Mx {TZx,0), because extrema (and saddles) lie 
at intersections of TZx and TZy. Accordingly, in normal 
speckle Mx {Tlx, 0) ~ Mx (7^, 0), Eqs. (|53|. In speckled 



speckle, however, the situation is very different. 

For vanishingly small K the only extrema are those of 
the a field, and the nodal domains of TZ are determined 
by these extrema. As K increases, however, numerous 
very small amplitude extrema associated with the pres- 
ence of the h field (here colloquially "6 extrema" ) begin 
to appear. These b extrema decorate with equal proba- 
bility both positive and negative nodal domains of 7?., but 
because of their very small amplitudes they do not divide 
up these nodal domains. Thus, Mq — Mx {TZ, 0) is not 
significantly increased by the presence of the b extrema. 
Ml — Mx {TZx , 0) , however, is increased in proportion to 
the square root of the density of these extrema, because 
TZx (and TZy) goes to zero at an extremum no matter how 
small its amphtude. As a result, A/1,0 = Mi /Mo begins 
to grow with K. This growth continues until the ampli- 
tude of the b source field is such that the amplitudes of 
the b extrema become large enough to carve up the nodal 
domains of TZ into small regions. When this happens Mq 




FIG. 12: Zero crossing densities iV^^'^', Eq. (|54b|, and 



N. 



(GG) 



Eq. 



(55b I, for two Gaussians (GG) for p = 0.01. 

2*1'^') is enhanced at its peak, 
10~* {K = 10~^), by a factor of 50 (44) over its normal 



For this value of p, iVi''^°^ (iV, 
K 



speckle value, Eqs. (53a 
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begins to grow rapidly. Thus the peak in Afifl is due 
to an initial growth of A/i with no significant growth in 
A/q, followed by a leveling ofi^ of A/i together with rapid 
growth in Mq. This scenario, which applies to all GDR 



source combinations, is illustrated in Fig. 13 



A/: 



2,1 



Why is A2 = Afx {T^xx,0) for speckled speckle so much 
larger than JVi — Afx {Tlx , 0), when for normal speckle the 
two are nearly equal, Eqs. ([SS])? The answer is analogous 
to the reason, discussed above, why in speckled speckle 
J\fi is so much larger than Ao, even though in normal 
speckle also these two densities are nearly equal. 

As discussed at the beginning of this section, for very 
small K each extremum is decorated on average by a 
single umbilic point, but as K increases the number 
of umbilics/extremum increases rapidly, Eqs. (47) and 

and 



Fig 



11 



Now, at an umbilic point TZxx — Tlyy 
[20]. Accordingly, as the number of umbilics 
increases so does J\fx {"R-xy, 0), the number of zero cross- 
ings of TZxy But for the circular Gaussian statistics 



xy 



relevant here N2 = ■\/5/3A/'x {TZxy,^), Eq. (50b I, so that 
as the number of umbilics increases so also does A2 . The 
net result is that because M2 increases rapidly with the 
number of umbilics, whereas A/i increases more slowly. 



only with the number of extrema (Fig. 13 1, A2/A/1 peaks 
at some intermediate value of K. These considerations 
could be illustrated by the analog of Fig. [13] in which 
A/q is replaced by A/i, M\ by A/2, and A/1,0 by N2,i] we 
refrain from presenting the resulting figure, however, be- 
cause except for scaling factors it is very similar to Fig. 



V. SUMMARY 

A priori, one might reasonably expect that adding a 
weak random field, the h field, onto a strong random field, 
the a field, would result in only minor perturbations of 
the a field structure and statistics. The opposite was 
found to be true, however, and just about every struc- 
tural and statistical property of the combined field was 
shown to display large anomalies when the optical power 
Pb of the b field is some lO"** - 10~^ the power Pa of the a 
field; these anomalies of speckled speckle are summarized 
below for p = a/b= 10"^ - 10"^. 

(i) Singularity clustering: Vortices in scalar fields and 
C points in vector fields cluster in the dark regions sur- 
rounding the bright speckle spots of the a field; the reason 
is that singularity formation requires equal amplitudes 
for the a and h fields. 

(ii) Phase field statistics: In both scalar and vector 
fields extrema can outnumber singularities by orders of 
magnitude; in contrast, in normal speckle singularities 
always outnumber extrema. 
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FIG. 13: (a) Zero crossings densities b ^Ni (thick curve) and 
b'^No (thin curve) vs. K for a and b fields Gaus sians (GG) 
with p = a/b = 0.01. A^o (Ni) is given in Eq. ( |49a| l, (Eq. 
(49b 1). The small white circles are the square root of the 
scaled density £ of extrema calculated from the exact r esul t 

y/£ (K) = (0F/3^/*) iVi (K), where £ is given in Eq. (44 1. 

The dotted curve is A''i,o = Ni/No rescaled to permit easy 
comparison with curves b~^No and b^^Ni. (b)-(d) Gray 
scale coded maps of TZ {x, y) for p = 0.01 for different values 
of K. In all maps TZ increases black to white, (b) K = 
lO""*. The thick curves are zero crossings of TZ that divide the 
field into positive/negative nodal domains. As can be seen, 
maxima (minima) almost invariably lie in positive (negative) 
domains. Maps for K < 10~* differ negligibly from the one 
shown here. (c) The central region in (b) enlarged by a 
factor of 100, so that the area shown is 10~* the area in (b). 
Thick white (black) curves are zero crossings of derivatives TZx 
{TZy). Every intersection of TZx and 7Zy marks a stationary 
point - extremum or saddle point; on average there are equal 
numbers of each type. Here the amplitudes of the extrema 
are so small that they are not apparent visually. (d) The 
region in (c) for K = 1. Thick curves are zero crossings of 
TZ. Here the extrema have large amplitudes. Because of 
this they are not only visible, but they also divide each nodal 
domain in (b) into some 10* small domains. Maps for K > 1 
differ unimportantly from the one shown here. 



(a) Scalar fields. For a sufficiently weak b field, say 
Pb/Pa ~ lO^'* - 10^^°, extrema and vortices maintain 
their normal ratio. But as Pb grows, the number density 
of extrema, E, grows faster than the density, V, of vor- 
tices, and the ratio E/V can become orders of magnitude 
larger than normal. When Pb/Pa ^ 1, the E/V ratio 
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returns to its normal value. Thus, E /V is normal for 
both very small and very large Ph/ Pa, and peaks at some 
still small, intermediate value of Pb/Pa- 

(&) Vector fields. For all nonzero Pf, less than some 
critical value P^, the ratio of the number density of az- 
imuthal extrema, E, can exceed the number density, V, of 
C points (Stokes vortices) by orders of magnitude. For 
Pb < Pbj E/V is essentially independent of Pb and is 
substantially equal to half the maximum scalar field ra- 
tio E/V for the same value of p; P^ approximately equals 
Pb at the peak of E/V. When Pb/Pa > 1, E/V returns 
to its normal value. Thus, E/V is large for small Pb and 
small for large Pb. 

(iii) Phase field derivative zero crossings: The ratio of 
the density of x-derivative zero crossings to y-derivative 
zero crossings changes anomalously with Pb', this anomaly 
is explained by the anomalous E/V ratio and the differ- 
ences in the geometry of the zero crossings at vortices 
and at extrema. 

(iv) Real and imaginary parts of the wave function: 

(a) Umbilic points. For a sufficiently weak b field, 
say Pb/Pa ~ 10"^" - 10^^^, extrema and umbilic points 
maintain their normal ^1:1 ratio. But as Pb is in- 
creased, the number density of extrema, £, grows faster 
than the density, U, of umbilic points, and the ratio 
£ /U can become orders of magnitude larger than nor- 
mal. £ /U returns to its normal value when Pb/Pa ^ 1- 

(b) Nodal dom,ains and derivative zero crossings. 
The densities of the zero crossings that define nodal 
domains and the densities of the zero crossings of all 
first and second order derivatives show anomalously large 
ratios. These anomalous ratios are explained by the 
anomalous £ /U ratio and the relationships of the differ- 
ent zero crossings to extrema and umbilic points. 

In general, the anomalous ratios of critical point den- 
sities occur because the number density of one wavefield 
feature grows rapidly with an increase in the power of 
the b field while the growth rate of a second feature lags 
behind. Such an imbalance in growth rates can be ex- 
pected to be the norm, and other properties of speckled 
speckle, including properties that defy calculation, can be 
expected to show large anomalies; these anomalies can, 
and should be, sought experimentally. 
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APPENDIX A: COMPOSITE FIELD 
AUTOCORRELATION FUNCTION 



Returning to the derivation of the autocorrelation 
functions of the simple sources in Eqs. ( 6a l-( 6c I we write 



D (p) of Eq. ^ as 

W (pr) = N (pr) /D (p) , (Al) 

where for the Gaussian (G), disk (D), and ring (R): 

7V(°) {pr) = 2p^exp {-p'^r^) , D^^' (p) = 2p^; (A2a) 
Ar(D) {pr) = p^J, (pr) / (pr) , D^^l (p) = p' /2- (A2b) 
7V(R) {pr) = epJo {pr) , {p) = ep. (A2c) 



Inserting Eqs. (A2al-(A2c) into Eq. (All immediately 
leads to Eqs. (|9a|)-(|9cp. 



The numerator (denominator) of the composite auto- 
correlation function in Eq. ([4| i s a linear combination of 
the N {D) in Eqs. ( |A2a )-( A2c I with p replaced by a or 
b, as required. Inserting the appropriate combinations of 
N and D into Eq. ^ leads to Eq. ([l2| with 



^(T„T,) ^ ^)/pf . 



(A3) 



with Pi"^"-* {P^^''^) the total optical power in source a 
(6). For a Gaussian (G), disk (D), and ring (R), 



p(D) _ 2t 

Pp(«) = 27rpe/p, 



(A4a) 
(A4b) 
(A4c) 



where p — a,b. For a Gaussian Ip is the peak intensity at 
the center of the source profile, for a disk Ip is the uniform 
source intensity within the disk, and for a ring Ip is the 
uniform source intensity within the narrow annulus. 

As a specific example, for a composite source with o a 
Gaussian and b a ring. 



exp (— r^)+K('^^)Jo{br) 



1 + if(GR) 
= belb/2a^la. 

APPENDIX B: DERIVATIVES 

The autocorrelation function W of H (or equivalently 
of J) is 

W{pr)^{n{0)n{pr))/{n^), (Bla) 
= {J {0)1 {pr)) /{!'), (Bib) 
= (7^2) W {pr) = (J^) W {pr) . (Blc) 

Suppressing momentarily for notational simplicity super- 
scripts T, 



r2n^{nl)/{n^), 

'd'^"W{pr) 



W in terms of the numerator N {pr) and denominator 



(B2a) 
(B2b) 



r=0 
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Reintroducing superscripts T, 



with 



1 + if(T„T,) 



where 



„(T) 



For the Gaussian (G), disk (D), and ring (R): 



„(G) 
-2 



2, cf ) 



12, ) 



120; 



1/4, C4 



(D) 
4 

(R) 



1/8, = 5/64; 

(R) 



1/2, cr = 3/8, 4"'^ 5/16 



(B3) 
(B4) 

(B5a) 
(B5b) 
(B5c) 



,(T) _ oJT) n) _ 



2c\'\ m'/> = (8/3)4^ \ 



(T) _(T) _ 



(16/5) 



,(T) 
'6 • 

(C5) 

4P is given in Eqs. |b5|) for a Gaussian (T = G), disk 
(T = D), and ring (T = R). 

As a specific example, for a composite source with a a 
Gaussian and b a disk, 

Mf °) = (384a6 + K(G^)byi) / (l + i^(GD)), 

Like the autocorrelation functions, the moments of 
composite sources are the optical power weighted aver- 
ages of the moments of the individual sources a and b. 



Defining derivatives for the annulus in Eq. (10a 



using Eq. ( B2b I , to leading order in the small quantity 
e/p, 

rfV4^^-l + (l/4)(£Vp'), (B6a) 

rfV4^^ = l + (5/6)(£Vp'), (B6b) 

4^)/4^) = 1 + (7/4)(£Vp^). (B6c) 

For e/p = 0.1, for the various quantities calculated in the 
text the differences between the experimentally realizable 
annulus and the theoretical ring are therefore less than a 
few percent, whereas for e/p — 0.01 these differences are 
of order 10^^. 

As a specific example for, say, r^, for a composite 
source with a a ring and b a disk, 

rf °) = {Says + K^^^hys) / (l + . 



APPENDIX C: MOMENTS 

The moments of the source function S are defined 
for a composite source by 

Expanding both W {pr) and Jq (ur) in Eq. 



{2n) 



I ' 2n 



It is useful to express M^f^''^ for speckled speckle in a 
form that parallels Eq. (Il2|, 



M 



where 



(T^Tfc) 
2n 



MZ^ya) + Ki™MZ;\b) 



1 + i^(T„Tt) 



M^P ip) = mg^p'-, 



(C3) 



(C4) 



APPENDIX D: 

We discuss here the dimensionless parameter rj that 
plays a central role in the level crossing statistics of 
derivatives of the phase, of the real and imaginary parts 
of the wavefunction, and of the intensity. Extending 
[17], we write for speckled speckle 



V 



:t„T,) _ ^(T„T,) 



/(4-»-^))^ 



1, 



(Dl) 



where the normalized derivatives ^'2^""'"'''' are given in Ap- 
pendix B. Of special interest is the regime where jyC^"'^'") 
becomes anomalously large. We find that jyC^"'^'') 
reaches a maximum, r]ma^''\ when xC^aTb) _ ICmax'''\ 
where 



(T„) 

ir(T„Ti,) ^ C2 2 

'^max (Xb) " ' 

^2 

/max 



1 



44^^)4'^^) \P 

For example, for a a ring (R) and b a disk (D) 

,(R) 



(D2a) 
(D2b) 



r(RG) _ !:2_ 2 _ 2n 



(RG) ^ 
/max 



„(D) 



44^)4°) VpV V 



For normal speckle, for a Gaussian 77 = 2, for a disk 
77 = 1, and for a ring 77 = ^. In contrast, for speckled 

speckle for say p — 0.1, r^nmx'' = 75, TyiiS.^'' = 50, and 
„(RR) _ 25 
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